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I Abstract 

We constrain the spectrum of primordial curvature perturbations V(k) by using 
the new data on the Cosmic Microwave Background (CMB) from the Boomerang and 
MAXIMA experiments. Our study is based on slow-roll inflationary models, and we 
consider the possibility of a running spectral index. Specifically, we expand the power 
\ spectrum P(k) to second order in In A:, thus allowing the power spectrum to "bend" 

J> ■ in fc-space. We show that allowing the power spectrum to bend erases the ability of 

0^ , the present data to measure the tensor to scalar perturbation ratio. Moreover, if the 

primordial baryon density fit,/i 2 is as low as found from Big Bang nucleosynthesis (BBN), 
the data favor a negative bending of the power spectrum, corresponding to a bump-like 
feature in the power spectrum around a scale of k = 0.004 Mpc -1 . 
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1 1 Introduction 

CO 

Inflation is generally believed to set the stage for the evolution of the universe, in particular 



providing the initial conditions for structure formation and cosmic microwave background 
(CMB) anisotropics. From a given inflationary model one can calculate the power spectrum 
of primordial curvature perturbations, V(k), which is a function of the wavenumber k. This 
power spectrum can then be Taylor expanded about some wavenumber ko and truncated 
after a few terms [fill 



InV(k) = lnV(k ) + (n - 1) ln£ + \4rT 

fco 2 a Ink 



ln 2 A + ... (1) 

k fe 



where the first term corresponds to a scale invariant Harrison-Zel'dovich spectrum, the 
second is the power-law approximation, and the third term is the running of the spectral 
index, which will be our main concern below. In order to provide the nearly scale-invariant 
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perturbations, which seem to be observed, one probably has to consider a slow-roll (SR) 
model for the later stage of inflation. The properties of SR models are well-known, and to set 
notation we recollect the basic features (see Ref. 0] for a review and list of references). In SR 
models one demands that the first few derivatives of the inflaton potential should be small. 
Traditionally this is expressed with the 3 SR parameters (e, ?7, £ 2 ), which roughly correspond 
to the first, second and third derivative of the potential. With these SR parameters one can 
express the scalar spectral index, n{k) — 1 = dlnV /dink, and tensor spectral index, ny(A;), 
and their derivatives M. We will here adopt a slightly different notation, and instead of 
the set (e,r/, £ 2 ) we will use the 3 parameters (n, r, d\ n k), where n = din V /dlnk\k=k + 1 is 
the scalar spectral index at the scale feo, the parameter r is the tensor to scalar perturbation 
ratio at the quadrupole, and d\ n k = dn/dlnk\k=k - The reason is simply that these 3 variables 
are closer related to what is being observed. With these 3 parameters one automatically 
expresses pi |I| the tensor spectral index and its derivative in k$ as 

(2) 



r driT r 

tit = and 



r 

(n-l) + — 
v ' 6.8 



6.8 dink 6.8 

The factor 6.8 in the above equation is actually model-dependent and should be calculated 
for each given model, in particular for different Q\ ||, but for simplicity we use the fixed 
value 6.8. 

The anisotropics, which are integrals over all the wavenumbers, pick up the major con- 
tribution from k ~ IHq/2, and hence one finds [Q] 

Ci [n(fc)] nyMi/io)d lnk 

Ci [n(ko)] * \TJ ' {3) 

where Iq ~ 2ko/Ho, which means that for a running spectral index, <9[ n £. ^ 0, the power 
spectrum will be bent (up or down depending on the sign of di n k), besides the normal tilt 
which arises for n / 1. In many SR models is expected to be very small since it is 
second order in the small parameters, however, there are very interesting models where this 
need not be the case |], f7|, and d\ n k may assume values big enough to be observable (see 
e.g. Refs. § |). 

Let us briefly mention where in our parameter space the different inflationary models 
lie. Traditionally |1C] one divides the simplest inflationary models into the following groups: 



small field (a < — 1), large field (—1 < a < 1), and hybrid models (a > 1), where a is defined 
through the SR parameters (see Appendix A for notation) 

r n — 1 

rj = at where e = , r\ = h 0.15 r . (4) 

13.6 2 

If we for simplicity consider models where the third derivative of the potential is zero 




Figure 1: The various slow-roll models presented in (n,d\ n k) space. The dotted line is the 
conservative limit, e < 0.1. The two dashed lines are the two attractors. 

= 0), then the different inflationary models are placed in the (n, 9i n / c )-space as shown 
in Fig. 1, where we have plotted various curves for different a. Here we have allowed only 
the conservative limit e < 0.1 (showed by the dotted line), and assumed that we only have 
to expand the expressions to leading order. Naturally the graph moves up and down by 
inclusion of the third derivative of the potential, £ 2 ^ 0. By changing the pivot scale, ko, 
around which the power spectrum is expanded, one can also move the graph sideways (see 
discussion in section ||), and the n in Fig. 1 should therefore be thought of as the n where 
the expansion scale has been chosen near the center of the probed scales. 



It is interesting that there seems to be attractors in the SR parameter space [11|. The 



two attractors found in [11] can be expressed as r = and r = — 6.8(n— 1), and the behavior 



of these solutions in the plane (n, d^k) is shown in Fig. 1 (dashed lines). As one can clearly 



see, according to these results d\ n k should be positive (or zero). Again, including the third 
derivative of the potential will also allow for a negative <9i n fc. From Fig. 2 one sees how these 
two attractors (dashed lines) are on the borders between the different inflationary models. 




Figure 2: The various slow-roll models in (n,d\ n k) space. The hybrid models correspond to 
a > 1, large fields means — 1 < a < 1, and small fields means a < — 1. The dashed lines 
are the two attractors. 

We will in the present paper discuss the ability of the present day data to provide 
information on the 3 parameters, (n, r, d\ n k)- 

2 The data 

Anisotropics in the CMB were detected for the first time in 1992 by the COBE satellite 
p2|. Recently, however, our knowledge of the temperature perturbations has increased 



dramatically with the results from the two balloon-borne experiments Boomerang [13] and 



MAXIMA [14]. One of the main conclusions from these experiments is the confirmation of 
the position of the first acoustic peak at I ~ 200, which strongly indicates a flat universe, 
^tot = 1- This seems to be a confirmation of the inflationary paradigm, since fitot = 1 is 
a rather general prediction of the simplest inflationary models. Very interestingly, the data 
also suggest that the second acoustic peak in the power spectrum is much less pronounced 



than predicted in flat models with baryon density compatible with BBN [16, 17 1. This could 
be an indication of new physics and accordingly a large number of papers on this subject 
have appeared since the release of the Boomerang and MAXIMA data. One possibility is 
that the primordial spectrum of fluctuations produced by inflation is not described by a 
smooth power-law, but rather that it has bumps and wiggles [18, 19, 20]. In the next section 



we will discuss this possibility further in light of our numerical results. 



The data from Boomerang and MAXIMA was recently analyzed in Ref. [15|, where 
a search in (n, r) space was performed. In that work KBl the assumption was made that 
dy n k = 0, which means that for each set of (n, r) one must carefully adjust the third derivative 
of the potential to make d\ n k = 0. This is perfectly possible, but it is more natural to allow 
d\ n k to vary, as we will do below. The effect of running of the tensor spectral index is very 
small, but for consistency we include it as described in Eq. (||). 

3 Data analysis 

In order to to investigate how the new CMB data constrain the parameter space (n, r, d\ n k) 



we have performed a likelihood analysis of the data sets from COBE [12], Boomerang |l 
and MAXIMA Q. The likelihood function to be calculated is 

£oc4exp^-£ j, (5) 

where i refers to a specific data point and 9 is a vector of cosmological parameters for the 
given model 

6 = {tt m , n A , fl b , H ,t, Q, n s ,r, d\ nk , ■ ■ ■} . (6) 

In the present case we have calculated the likelihood function for the following parameter 
space: O m , the matter density, f2b, the baryon density, Hq, the Hubble parameter, as 
well as the inflationary parameter space n,r,d\ n k- We have assumed that the universe is 
flat Q\ = 1 — O m , and that reionization is not important, r = 0. We also treat the overall 
normalization, Q, as a free parameter. The experimental groups quote estimated calibration 
errors for the experiments. These should also be accounted for, and we do this by allowing 
the data points to be shifted up or down by this amount: 10% for Boomerang [Q3[ and 



4% for MAXIMA [14]. It should also be noted that the data points are to some extent 
correllated. The likelihood function, Eq. (|5|), is based on the assumption that errors are 
uncorrellated, and therefore a (quite small) error is introduced into our analysis. However, 
until the full data sets from the experiments are publicly available, the magnitude of the 
error is difficult to quantify. 

Finally, we have chosen the pivot scale in Eq. (|l|) as ko = 0.05 Mpc -1 . This choice 
is made for convenience, since ko = 0.05 Mpc -1 is the scale at which wavenumbers are 
normalized in the CMBFAST code. Our main results, however, do not depend on the 
specific value of ko. 

In Fig. 3 we show the allowed region in the inflationary parameter space for the case 
corresponding to the calculation of Ref. |jl5[ |, namely d\ n j~ = 0. In the left panel the allowed 
region was derived assuming that Qi,h 2 = 0.019 (BBN prior plf|), and the left panel of our 
Fig. 3 thus corresponds to their Fig. 3. Our result is almost identical to theirs, the only 
difference being that they have used a slightly larger space of other cosmological parameters, 
so their constraints are slightly less restrictive than ours. In Fig. 4 of Ref. [|l5fl , r was allowed 
to vary. Changing the value of r does change the likelihood contours (pushing them to larger 
values of n), but only to a small extent. This illustrates that allowing r as a free parameter 
in our analysis would lead to slightly larger allowed regions, but would not in any way 
invalidate our conclusions. 

One thing which should be noticed is that even the best fit point has x 2 /d.o.f. = 1.2, 
indicating that it is a quite poor fit. This is not too surprising because it is well known that 
no good fit to the new CMB data can be obtained with a low baryon density, even allowing 
r to vary (if], [l7| . 

In the right panel we show the allowed region for the case where f^/i 2 = 0.030, which is 
the value favored by the CMB measurements. In this case we find an allowed region very 



similar to that found in Fig. 2 of Ref. [15|, although again our allowed region is slightly 
smaller than theirs because we use a smaller parameter space. It should also be noticed 
that the best fit point now has x 2 /d.o.f. = 0.53, a very good fit. This corresponds well to 



the findings of other likelihood analyses [16, 17], that the CMB data can be very well fitted 
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Figure 3: The allowed region in the n, r parameter space, calculated from the combined 
COBE, Boomerang and MAXIMA data. The dark shaded (green) regions are la and the 
light shaded (yellow) are 2a. The left panel assumes a BBN prior on Q^h 2 = 0.019, whereas 
the right panel is for ^h? = 0.030, the value which best fits the CMB data. Note that the 
best-fit points, marked by diamonds, are really at r = 0, but have been shifted slightly so 
that they are more visible. 



in models with high baryon density. 

If the assumption d\ n k = is relaxed, the results change substantially. In Fig. 4 we 
show results for the likelihood analyses where d\ n k is allowed to vary. Again, the left panels 
correspond to Qbh 2 = 0.019 and the right panels to O^/i 2 = 0.030. As one can see, we have 
varied the inflationary parameters in a wide range, even larger that what allowed by SR 
approximation. The reason is that we believe that the main conclusions that we obtain are 
quite general and give indications on the shape of V(k) which are relevant even outside the 
context of SR models. 

In the upper panels of Fig. 4 we show the same allowed region as in Fig. 3, when the 
assumption d\ n k = is relaxed. When d\ n k is allowed to vary, the preferred region in 
parameter space is shifted completely. The best fits are still for a model with no tensor 
component (for both values of Sl^h 2 ), but now there is no real constraint on r. The best fit 
models now have x 2 /d.o.i. = 0.54 and 0.57 respectively, which indicate very good fits. 
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Figure 4: The allowed region in the n, r, d\ n k parameter space, calculated from the combined 
COBE, Boomerang and MAXIMA data. The dark shaded (green) regions are lcr and the 
light shaded (yellow) are 2a. The left panels assume a BBN prior on l^/i 2 = 0.019, whereas 
the right panels are for f^/i 2 = 0.030, the value which best fits the CMB data. Note that the 
best-fit points, marked by diamonds, are really at r = 0, but have been shifted slightly so 
that they are more visible. The solid line in the left (n, d\ n k) panel corresponds to the power 
spectrum exhibiting a distinct feature at k = 0.004 Mpc -1 (see Eq. (8) and the surrounding 
discussion) . 
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Figure 5: Power spectra for four different models. Panel (a) is the best fit model with 
dink = 0, r = 0, (b) the best fit with d\ n k = 0, r = 2, (c) the best fit with d\ n k ^ 0, r = 0, and 
(d) the best fit with d\ n k ^ 0, r = 2. The data points are from the Boomerang experiment 
flUI . The curves show: The scalar component (dotted lines), tensor component (dashed 
lines), and the combined fluctuation spectrum (solid lines). 

In order to understand why allowing d\ n k to vary erases any ability to constrain r, it is 
instructive to look at the power spectra for some of the best fit models directly. For the 
remainder of this section we will assume that fl^h 2 = 0.019, in accordance with Big Bang 
nucleosynthesis. In Fig. 5 we show four different power spectra, all calculated for the case of 
Qbh 2 = 0.019. Panels (a) and (b) are both for d\ n k = 0. Model (a) is the best fit with r = 0, 
whereas (b) is the best fit with r = 2. Model (b) is a very poor fit because introducing a 
tensor component while still fitting the high /-values severely underestimates power around 
the first peak. This cannot be remedied by shifting n alone. Therefore, the allowed values 
of r are quite tightly constrained. 

Panels (c) and (d), on the other hand, show models where d\ n k is allowed to vary. Model 
(c) is the best fit for r = and (d) is the best fit for r = 2. In contrast to the case with 
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dink = 0, it is possible to obtain a decent fit, even with a large tensor component because 
the power spectrum can be "bent" by having a non-zero d\ n k- 

2.5 r i 




-4 -3 -2 -1 

10 10 10 10 1 



k - Mpc " 1 

Figure 6: Power spectra of scalar perturbations for selected values of the parameters (n, d\ n k) 
which provide good fits to the CMB data. The solid line corresponds to the best-fit point, 
(n,d\ n k) = (0.5, —0.2). The dashed lines correspond, from top to bottom in the maximum 
region, to (n,d lnk ) equal to (0.4,-0.24), (0.6,-0.16), (0.7,-0.12), respectively. The dot- 
dashed line corresponds to best-fit model with constant spectral index, i.e. (n,c\ n k) = 
(0.9,0). All normalizations are arbitrary. 

The effect of a running spectral index can also be understood by looking at Fig. 6, where 
we show the power spectra of scalar perturbations corresponding to values of (n, <9i n fc) which 
provide good fits to the CMB data. The effects of a d\ n k ^ is somewhat similar to introduce 
a feature in the power spectrum at a fixed scale. Specifically, for d\ a k < 0, the power spectra 
show a maximum at a scale k m given by [] 

m(& m /A!o) = -(n - l)/b\ nk . (8) 

The CMB data provide a constraint for the possible positions of this maximum. This can 

be understood by looking at the left (n, c\ n k) panel of Fig. 4, from which it is evident that 

the allowed region in the plane (n, <9i n fc) lies along a line corresponding to k m = const. 
4 One can easily show that the power spectrum (p.h can be rewritten as a Gaussian in In k, centered around 

' — ■ 1/2 

km = ko exp [— (n — l)/9i n fe] and having a width equal to d lnk , i.e.: 

InP(fe) =lnP(fc m )+ ^ ln 2 -^ + --- (7) 
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Specifically, the best fit model corresponds to k m = 0.004 Mpc , while possible values are 
k m = 0.0015 - 0.01 Mpc" 1 . 

It is easy to understand why such a feature helps to fit the CMB data. In order to 
increase sizeably the first to second peak ratio in the CMB, one needs a power spectrum 
with a large negative slope at intermediate and small scales, say e.g. k > 0.01 Mpc -1 . If 
we assume d\ n k = this is not allowed; the power spectrum is then a monotonic function 
of k and, as consequence, one automatically obtains an excess of power at large scales, 
corresponding to COBE normalization. If instead one has d\ n k < 0, the power spectrum is 
non-monotonic and has a maximum at a specific scale k m . Specifically, if k m ~ 0.004 Mpc -1 , 
the power spectrum decreases both at large and at small scales so that one has no trouble 
in reproducing the COBE data. Clearly, as a by-product, one looses the ability to constrain 
the tensor to scalar ratio r. The scalar power spectrum at large scale can in fact be strongly 
suppressed and this allows for a large contribution from tensor perturbations. 

We conclude this section by comparing our result with the result of [jig] , in which a 
Gaussian bump in logfe was added to a standard power spectrum and its position k\, was 
constrained by CMB data. In very nice agreement with our results they found a best fit fc& 
of roughly 0.005/i Mpc -1 , and an allowed region for ky, of O.OOl/i — O.Ol/i Mpc -1 . 

4 Conclusion 

We have considered the ability of the present day CMBR data to distinguish between various 
slow-roll models of inflation, allowing the scalar spectral index, n{k), to vary with scale. 
Specifically, we have expanded the power spectrum V(k) to second order in In A: and we 
have derived the constraints on the parameter space (n, r, d\ n k) which can be obtained by 
using COBE, Boomerang and MAXIMA data. 
We have seen that: 

i) If we allow d\ n k ^ the tensor to scalar ratio r is essentially unconstrained, even if we 
assume Clbh 2 as low as suggested by BBN considerations. 

Moreover, assuming a BBN prior of ^h 2 = 0.019, we have found that: 

ii) A negative bend of the power spectrum, d\ n k < 0, is favoured by the CMB data. 
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iii) The best fit model, (re, dink) = (0-5, —0.2), which provides a very good fit to the CMB, 
corresponds to a power spectrum which deviates quite strongly from a power law approxi- 
mation. The large values obtained for re — 1 and d\ n k are bordering to invalidate the slow-roll 
approximation, however, one should keep in mind that changing the pivot scale, ko, will 
change the value of re — 1, and inclusion of the 3 rd derivative of the inflaton potential will 
change the value of d\ n k; 

iv) The CMB data favor power spectra with a bump-like feature at scales k m = 0.0015 — 
0.01 Mpc~\ the best fit value being k m = 0.004 Mpc -1 . 

Summarizing, the general result of our analysis is that a single power law for V{k) does 
not provide a good fit to CMB data, if we assume a BBN prior f^/i 2 = 0.019. In particular, 
CMB data favor models in which V(k) deviates strongly from a power law approximation. 
Specifically, the scale at which the power law approximation should be broken, i.e. second 
(or higher) order terms in eq. ([!]) become important, is k = 0.0015 — 0.01 Mpc -1 . 
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A Notation 



We use the notation: 

M 2 fV'\ 2 2 V" M 2 fV'\ 2 

see |jl5| for details. The notation with r\ = M 2 V" /V, used e.g. in [Q], simply corresponds 
to the substitution a — > a + 1 in eq.@. Independently on the chosen definition for rj, one 
finds 



9i„ fc = -2e 2 + 4— 

0.0 



(n - 1) H 

v y 2 6.8 



(9) 



where £ 2 = M V / V'"/V 2 . One could potentially include higher order terms in Eq. (1) 
corresponding to d 2 n/dhik 2 ^ 0, which would be expressed as 

d 2 n S „ 3 , 1 * ( n r 



**> ^ + -^ 9- -(„-!)] (10) 



^ 2 / \ 3 



-2(n-ir 15(n-l) — -15 — , (11) 

V ; 6.8 V ' V6.87 V6.87 V ; 

|^ = d lnk ^--(n-lf^--3(n-l)(^-) 2 -2(^-) 3 , (12) 
dlnA; 2 6.8 V ; 6.8 V \6.8/ V6.8y ' V ; 

with d 3 = 2eV""/V. 
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